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We report the first search for CP violating decays of the T(4S) using a data sample that contains 
535 million T(4S) mesons with the Belle detector at the KEKB asymmetric-energy e^e~ collider. A 
partial reconstruction technique is employed to enhance the signal sensitivity. No significant signals 
were observed. We obtain an upper limit of 4 x lO^'^ at the 90 % confidence level for the branching 
fractions of the CP violating modes, T{4S) B°B° J/i^K^ + J/tp{ric)K^- Extrapolating the 
result, we find that an observation with 5a significance is expected with a 30 ab~^ data sample, 
which is within the reach of a future super B factory. 
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CP violation has been established in the neutral kaon 
system [l| and the neutral B meson system In 
the standard model (SM) Kobayashi-Maskawa theory, it 
arises from an irreducible phase in the weak interaction 
quark-mixing matrix [1]. This theory predicts that CP 
violation in the T(45') system should also exist. 

In the decay T(4S') ^ B°B° ^ /1/2, where /i and 
/2 are CP eigenstates, the CP eigenvalue of the final 
state /1/2 is ^ = —^i£,2- Here the minus sign corresponds 
to odd parity from the angular momentum between fi 
and f2- If /i and /2 have the same CP eigenvalue, i.e. 
{^1,^2) = (+1,+1) or (-1,-1), equal to -1. Such 
decays, for example (/i,/2) = {J/i^Kg, J/^jjKg), violate 
CP conservation since the T(45') meson has J^^ = 1 
and thus has ^t(4S) = +1- The branching fraction within 
the SM is 

B{ T{AS) ^ ^ /1/2 ) 

= F ■ B{T{4S) B°B°)B{B" fi)B{B° /2), (1) 

where _F is a suppression factor due to CP violation. The 
factor F can be calculated in terms of mixing and CP 
violating parameters [3], 

F ^ ^1^(2 sin 201 f (2) 
1 + 

= 0.68 ±0.05, 

where x = Amd/T = 0.776 ± 0.008 @, Arud is the B° 
mixing parameter, F is the average decay width of the 



neutral B meson. The angle (fii is one of the three inte- 
rior angles of the unitarity triangle of the quark-mixing 
matrix, and sin2(/)i = 0.675 ± 0.026 ^. The effect of 
direct CP violation is neglected in this formula. The 
same expression also holds for the case in which /i and 
/2 are different final states both of which are governed 
by 6 — > CCS transitions; examples include rjcKg, Tp(2S)Kg 
and XciA'g. 

In this Letter, we present the first search for CP vio- 
lating decays of the T(4S'). The data sample used con- 
tains 535 million T(4S') mesons collected with the Belle 
detector at the KEKB asymmetric-energy e~^e~ (3.5 on 
8 GeV) collider 0. The Belle detector is a large-solid- 
angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber 
(CDC), an array of aerogel threshold Cherenkov counters 
(ACC), a barrel- like arrangement of time-of-flight scintil- 
lation counters (TOF), and an electromagnetic calorime- 
ter comprised of CsI(Tl) crystals (ECL) located inside a 
superconducting solenoid coil that provides a 1.5 T mag- 
netic field. An iron fiux-return located outside of the 
coil is instrumented to detect mesons and to identify 
muons (KLM). The detector is described in detail else- 
where [7|. Two inner detector configurations were used. 
A 2.0 cm radius beampipe and a 3-layer silicon vertex 
detector were used for the first sample of 152 million BB 
pairs, while a 1.5 cm radius beampipe, a 4-layer silicon 
detector and a small-cell inner drift chamber were used 
to record the remaining 383 million BB pairs 
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The identity of each charged track is determined by a 
sequence of hkehhood ratios that determine the hypoth- 
esis that best matches the available information. Tracks 
are identified as pions or kaons based on their specific 
ionization in the CDC as well as the TOF and ACC re- 
sponses. This classification is superseded if the track is 
identified as a lepton: electrons are identified by the pres- 
ence of a matching ECL cluster with energy and trans- 
verse profile consistent with an electromagnetic shower; 
muons are identified by their range and transverse scat- 
tering in the KLM. 

We use 2.68 x 10^ Monte Carlo (MC) simulation events 
for each signal category. For background MC events, wc 
use a sample of 3.9 x 10^° generic BB decays in which one 
of the B mesons decays to a known J/ip(fi~^ii~ or e^e~)X 
final state. For the datasct used in the present anal- 
ysis, the MC simulation predicts a small signal yield, 
0.04 events, when we choose the combination (/i, /2) = 
[J/ipKl, J/i'K°) and fully reconstruct both J/V'A'g fi- 
nal states. Here we use the J/ip — > e+e~, fJ,~^n~ and 
Kg — > 7r+7r~ modes. In order to increase the signal 
yield, we instead adopt a partial reconstruction method. 
We fully reconstruct one 5° J/ipKg decay (called 
fj/i>K° hereafter) and find another K^g (called '^^j^o 
hereafter) from the remaining particles. We then recon- 
struct the recoil mass (Af'^'^"") using J/ii)Kl and '^^SA^g. 
The recoil mass distribution should in principle include 
peaks that correspond to the rjc-, J/'iA, Xcii or -0(25'). 
We choose two of the possible combinations, (/i, /2) 



[fj/4,K0, J/^^^'^Kfj and (^fj/^Kp rj^ 



tag^O 



In 



the following, these are referred to as inclusive- J/i/" com- 
binations and an inclusive- t/c combinations, respectively. 
Based on a MC study, we expect that the signal yield will 
increase by a factor of 40 compared to full reconstruction 
while maintaining a reasonable signal to background ra- 
tio (S/B) of about 1/7 for these two combinations. We 
do not use other combinations because the S/B ratio is 
less than 1/100. 

We use oppositely charged track pairs to reconstruct 
J/'f/' — > e+e", fJ.^ ^J.' decays, where at least one track is 
positively identified as a lepton. Photons within 50 mrad 
of the e"*" and e~ tracks are included in the invariant 
mass calculation (denoted as e~^e~{j)). The invariant 



mass is required to lie in the range 
M„e(^) - TOj/^ < 0.036 GeV/c^ and 



I- J ftp 



-0.15 GeV/c2 < 
-0.06 GeV/c^ < 
^j/ip < 0.036 GeV/c^, where toj/^ denotes the 
nominal mass of J/^p, Mee{j) and M^^ are the recon- 
structed invariant masses from e^e~['~f) and fi'^ fi~ , re- 
spectively. Asymmetric intervals are used to include part 



of the radiative tails. 



Candidate Kg 



decays 



are oppositely charged track pairs that have an invari- 
ant mass within ±0.016 GeV/c^ (~ 4a) of the nominal 
mass. The tt+tt^ vertex is required to be displaced 
from the interaction point in the direction of the pion 
pair momentum for *^^Kg. 
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FIG. 1: Mbc (left) and AE (right) distributions for B° 
J/ip{l+r)K°g{Tv+TT-) decay {I = e,fj,). The solid curves 
show the fits to signal plus background distributions, and the 
dashed curves show the background distributions. 



For the full reconstruction of a i? decay, we use the 
energy difference AE' 



TTicms 



-^b™am ^^"i beam- 



energy constrained mass il/bc 



l^bcam) - (Pb ) 



where i?b"am beam energy in the center-of-mass 

system (cms) of the T(4S') resonance, and E'j^^ and 
p^** are the cms energy and momentum of the recon- 
structed B candidate, respectively. The Mbc and AE 
distributions are shown in Fig. [T] The signal is extracted 
from an unbinned extended maximum likelihood fit to 
the Afbc-Ai? distribution. The signal shape is modeled 
with a single (double) Gaussian while the background 
shape is modeled with an ARGUS function [§] (a first 
order polynomial) for the Mbc {AE) distribution. We 
obtain 8283 ± 94 f j /^k° events when we do not require 
a ^""^Kl. 

We require 5.27 GeV/c^ < Mbc < 5.29 GeV/c^ and 
\AE\ < 0.04 GeV for fj/^K%- The recoil mass is calcu- 
lated by combining a fj/^K° candidate and a ^^^K^ can- 
didate. The expected number of signal events estimated 
from MC is 1.1 (0.6) with a reconstruction efficiency 
of 28.8 (26.8) % for the inclusive- J/i/; [rjc) combination 
where branching fractions of sub-decays are not included. 
With the partial reconstruction technique, the number of 
J/'(/' e+e", fJ,~^fJ,~ decays in the {J/ipKg, J/tpKg) 
combination is about twice as large as that for the 
{J/ipKg, rjcKg) combination. A total of 1.7 signal events 
are then expected in our datasct. 

The dominant source of background is generic B'^ de- 
cays. A partially reconstructed B candidate should be 
flavor non-speciflc if it is a signal event. On the other 
hand, about a half of the generic 5° decays that sur- 
vive the selection are flavor specific. In order to distin- 
guish between the signal and the background, we there- 
fore identify the flavor of the partially-reconstructed ac- 
companying B meson using leptons, charged pions and 
kaons that are not associated with the fully reconstructed 
B meson. The procedure for flavor tagging is described 
in Rcf . [l3| • We use an event-by-event flavor-tagging di- 
lution factor, r, which ranges from r = for no flavor 
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discrimination to r = 1 for perfect flavor assignment. 

We determine the signal yield by performing an un- 
binned extended maximum-likelihood fit to the candidate 
events. The likelihood function is 



^=^cxp 



N 



E-^- n 



E 



^ -f / A rl'ccoil _ 



1 L fe 



(3) 



where is the total number of candidate events, is the 
number of events and fk is the probability density func- 
tion (PDF) for each event category fc, which is inclusive- 
J/V", inclusive-77c or background. The parameters M™'^°'' 
and ri are the recoil mass and r value for the i-th event. 
The PDFs are obtained from the MC simulation. The 
recoil mass distributions are modeled with a triple Gaus- 
sian for each signal mode and an exponential shape for 
background. We do not find any peaking background in 
either the MC samples or in the M\,c sideband data. The 
PDFs for the r distributions are histograms with 10 bins 
obtained from MC. The ratio between the inclusive- J/ t/i 
and Tjc signals is fixed from the MC. 

We check the method using charged B de- 
cay control samples, T(4S') B^B^ 
{fB+, Jli> ^"^^K^ and r]c ^^^^K'), where /b+ stands for 
J/ip{e+er,^+^')K+ and D°{K+tt~ ,K+TT-iT+TT-)Tr+ 
decays Figure [2] shows the recoil mass distribution 
for the charged B control samples. The fit yields 206 ±57 
signal events, which is in good agreement with the MC 
expectation (183 events). If we float the ratio between 
the inclusive- J/V' and rjc modes, we obtain 96 ± 23 and 
109 ± 25 events for the inclusive- J/^/i and rjc modes, 
respectively. These results are also consistent with the 
MC expectation, 90 (93) events for inclusive- J/V' [rjc) 
mode. We obtain correction factors, the mean and width 
for the signal peaks and the slope for background, by 
fitting these samples. 

We adopted a blind analysis method and estimated 
systematic uncertainties before obtaining the final result. 
The systematic uncertainties for the combined branching 
fraction, B {T{AS) B°B° J/ipK°, {J/iP, r]c)K°) , 
are summarized in Table ID The dominant source of 
systematics is due to the uncertainties in the correction 
factors for the recoil mass distribution; we assign 20.5 %, 
which is the sum in quadrature of 19.7 % from the signal 
shapes and 5.5 % from the background shape. 



TABLE I: Systematic uncertainties in the branching fraction 
measurement. 



Source 


(%) 


Recoil mass distribution 


20.5 


r distribution 


4.2 


Reconstruction efficiency 


5.7 


Number of BB pairs 


1.3 


Branching fractions of sub-decays 


10.9 


Total 


24.3 
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FIG. 2: Recoil mass distribution for the charged B decay con- 
trol samples (top), recoil mass (middle) and r (bottom) dis- 
tribution for the neutral B decay control samples. The solid 
curve shows the fit to signal plus background distributions 
while the dashed curve shows the background distribution. 



Possible diff'erences between data and the MC in 
the r distributions are also studied. We use neu- 
tral B decay control samples, T(4S') B^B'^ 
(/bo I [J I4't11c) ^^^K^^ decays, where /^o represents 
B^ D^*'>^iT^ and D*~ followed by the decays 
D*- D°TT-, D° K+TT-, K+TT- 
/\+7r^7r^, /9+ tt+tt" and tt*^ 
35±16 signal events for these samples, which is consistent 
with the MC prediction (64 events) within two standard 
deviations. There is no discrepancy between data and 
fit results either in recoil mass or in the r distributions, 
as shown in Fig [2l We repeat the fit using the back- 
ground r PDF determined from the data in the recoil 



77. We obtain 



mass sideband regions Af^™'' G (2.40, 2.85) and (3.20, 
3.30) GcV/c^. The difference between the two fit results 
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FIG. 3: Recoil mass (upper) and r (lower) distribution for 
samples reconstructed as T{AS) (^J/ipKg , {J /tp,r)c)K's) 
decay. The solid lines show the fits to signal plus background 
distributions while the dashed lines show the background dis- 
tributions. 



at the 90 % confidence level, where the SM prediction is 
1.4 X 10"'^. This corresponds to < 2 at the 90 % 
confidence level. We also search for {j/tpKg, J/^Kg) 
combinations by fully reconstructing both B mesons. No 
candidates are observed. 

In summary, a search for CP violation in T(4S') de- 
cays was performed. In a data sample of 535 mil- 
lion BB pairs obtained via decays of the T{AS) reso- 
nance, no significant signals were observed. We obtain 
an upper limit of 4 x 10~^ at the 90 % confidence level 
for the branching fraction of the CP violating modes, 
T(4S') B"B" J/i^Kl + {Jjii, T]c)K%. Assuming the 
SM, with an integrated luminosity of 30 ab~^ that is ex- 
pected to be available in a future B factory, these decays 
can be observed with 5cr significance. 
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the Nil for valuable computing and Supcr-SINET net- 
work support. We acknowledge support from MEXT and 
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KIP of CAS (China); DST (India); MOEHRD, KOSEF 
and KRF (Korea); KBN (Poland); MES and RFAAE 
(Russia); ARRS (Slovenia); SNSF (Switzerland); NSC 
and MOE (Taiwan); and DOE (USA). 



(2.6 %) is included in the systematic error from the r 
distribution. We also repeat the fit without using the r 
distribution, which yields a result that differs by 3.3 % 
from the nominal fit result. We assign a 4.2 % system- 
atic uncertainty for the r distribution, which is the sum 
in quadrature of these two errors. 

Systematic uncertainties from event reconstruction are 
studied by varying the particle identification, A'^ selec- 
tion and other requirements. The resulting changes in 
the signal yield in data and MC for B^ J/ijiK^g and 
_B+ J/ipK^ are used to estimate the systematic error. 
In total, 5.7 % of the systematic uncertainty that is ob- 
tained from the sum in quadrature of differences between 
data and MC is assigned for event reconstruction. The 
uncertainty in the total number of BB pairs is 1.3 %. 
Uncertainties in the daughter branching fractions (H| are 
dominated by those for the 77c decays. 

The results of the final fit are shown in Fig. [31 
The extracted signal yield, —1.5 I2 g events, is con- 
sistent with zero as well as with the SM predic- 
tion (1.7 events). An upper limit is determined 
with a frequentist method [I4], where the PDFs are 
smeared to include systematic uncertainties. We obtain 
BiTiAS) ^ B°B° ^ J/i:K°,{J/i:,r]c)K°) < 4x10"^ 
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